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Abstract We consider stochastic differential systems driven by continuous 
semimartingales and governed by non-commuting vector fields. We prove that 
the logarithm of the flowmap is an exponential Lie series. This relies on a 
natural change of basis to vector fields for the associated quadratic covariation 
processes, analogous to Stratonovich corrections. The flowmap can then be 
expanded as a series in compositional powers of vector fields and the logaritm 
of the flowmap can thus be expanded in the Lie algebra of vector fields. Further, 
we give a direct self-contained proof of the corresponding Chen-Strichartz 
formula which provides an explicit formula for the Lie series coefficients. Such 
exponential Lie series are important in the development of strong Lie group 
integration schemes that ensure approximate solutions themselves lie in any 
homogeneous manifold on which the solution evolves. 
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1 Introduction 

We are concerned with Ito stochastic differential systems driven by continuous 
semimartingales and governed by non-commuting vector fields of the following 
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form 

d r>t 

Y t=Y 0 +y\ Vi(Y r ) dX' T , 

i=i J o 

for time t e [0,T] for some T > 0. Here the solution process Y t is K^-valued 
for some TV £ N. For each i = 1 ,,d, the X 7 t are driving scalar continuous 
semimartingales, and associated with each are governing vector fields Vi which 
are sufficiently smooth and in general non-commuting. Our goal herein is to 
compute the logarithm of the flowmap for such a system, i.e. the exponential 
series for the flowmap, and establish that it is a Lie series. The exponential 
series for the flowmap for stochastic differential systems driven by general 
continuous semimartingales was derived in Ebrahimi-Fard, Malham, Patras 
and Wiese HJ. What we achieve that is new in this paper is we: 

1. Establish the abstract algebraic structures that underly the flowmap and 
computation of functions of the flowmap in the context of general contin¬ 
uous semimartingales; 

2. Show by a suitable change of coordinates, the exponential series is a Lie 
series; 

3. Give a direct self-contained proof of the corresponding Chen-Strichartz 
formula which provides an explicit formula for the Lie series coefficients. 

The key idea that underlies establishing the exponential series as a Lie series is 
to express the flowmap in Fisk Stratonovich form for which the standard rules 
of calculus apply; see Protter |33 . A crucial integral ingredient in this step is 
that the Fisk-Stratonoviclr formulation of the flowmap can be expanded in a 
basis of terms involving solely compositions of vector fields—without any sec¬ 
ond order partial differential operators. We can then compute the logarithm 
of the Fisk -Stratonovich representation of the flowmap. This can be accom¬ 
plished in principle via the classical Chen -Strichartz formula using the shuffle 
relations satisfied by multiple Fisk-Stratonovich integrals as well as utilizing 
the Dynkin-Friedrichs-Specht-Wever Theorem to expand the logarithm in Lie 
polynomials of the vector fields. We subsequently convert the multiple Fisk- 
Stratonovich integrals back into multiple Ito integrals. This procedure thus 
generates an Ito exponential Lie series. That the logarithm of the flowmap 
is in fact an exponential Lie series is important for example, for the devel¬ 
opment of strong stochastic Lie group integration methods. See Malham and 
Wiese m for the development of such methods for Stratonovich stochastic 
differential systems driven by Wiener processes, for example those based on 
the Castcll-Gaines numerical simulation approach, see Castcll and Gaines 0] 

The development of exponential solution series for deterministic systems 
originates with the work of Magnus [5S] and Chen [B] in the 1950’s, and more 
recently with Strichartz [30] • Its development and early application to stochas¬ 
tic systems is represented by the work of Azencott [ I], Ben Arous [2], Castell 
and Gaines 01 El and Baudoin [3. Also see Fliess 03] and Lyons [25] for its 
development in control and theory of rough paths, respectively. The shuffle 
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product was cemented in firm foundations by the work of Eilenberg and Mac 
Lane [12] and Schiitzenberger [35], also in the 1950’s. The quasi-shuffle prod¬ 
uct is a natural extension of the shuffle product, for example to multiple Ito 
integrals. For a selective insight into its recent development in this context, 
see Gaines mm, Hoffman 03, Ebrahimi Fard and Guo [9], Hoffman and 
Ihara [U; and Curry, Ebrahimi-Fard, Malham and Wiese 0. 

Our paper is structured as follows. In Section0we derive the Ito chain rule 
and flowmap for systems driven by continuous semimartingales. Then in Sec- 
tion[3]we establish the abstract algebraic structures that underpin the flowmap 
and its logarithm. We endeavour to keep the connection to the stochastic dif¬ 
ferential system of interest and provide illustrative examples. We define the 
transformation to Stratonovich form we require in Section [I] and prove that 
the exponential series is a Lie series. Our direct self-contained derivation of 
the Chen-Strichartz coefficients is provided in Section [5] Lastly, we provide 
some concluding remarks in Section [6] 


2 Ito chain rule and flowmap 

Consider the Ito stochastic differential system given in the introduction of the 
form 

Yt =Y 0 +y f Vi(Y r ) dX*. 

i=i 

Here for each i = 1,... ,d, the X[ are driving scalar continuous semimartin¬ 
gales on a filtered probability space (l7, T. P) satisfying the usual con¬ 

ditions of completeness and right-continuity. We assume without loss of gener¬ 
ality that the X 1 are chosen such that the quadratic covariations [X 1 , X 1 ] = 0 
for all i # j. The V) are associated governing vector fields which we assume 
are sufficiently smooth and in general non-commuting. We suppose the solu¬ 
tion process Y t) which is R^-valued for some N e N, exists on some finite 
or possibly infinite time interval. In coordinates the vector fields V. for each 
i = 1 ,...,d act as first order partial differential operators on any function 
/: 1^ —» M as follows 


N 

Vi : /00 - 2 Vi(Y)d Yj f(Y). 

3 = 1 

For brevity we will often express this vector field action as {Vi • d)f(Y) or 

VofoY. 

Definition 1 (Flowmap) We define the flowmap <p t as the map prescribing 
the transport of the initial data f oYq to the solution / o Y t at time t for any 
smooth function / on T8L N , i.e. ipt '■ f ° Yq h-» / o Y t . 

The solution Y t = ip t o id o F 0 corresponds to the choice / = id, the identity 
map. The Ito chain rule is the key to developing the Taylor series expansion 
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for the solution Y t about the initial data. The Ito chain rule implies that for 
any function /: 1^ —> M, the quantity f(Y t ) satisfies 

f(Y t ) = f(Y 0 ) + Y l \\v i -d)f(Y r ) (LY‘+ i 2 d 2 )f(Y T )d[X\X'] T , 

i= i J o i= i J o 

see for example Protter m • In this formula we have used the notation 

N 

{V®V.:d 2 )f(Y):= £ V>(Y)Vf(Y)d Yj d Y J(Y), 

j,k=l 


while for each i = 1,... ,d the terms [X l ,X l ] represent the quadratic vari¬ 
ation of X 1 . At this stage it makes sense to extend, first our set of driving 
continuous semimartingales to include these quadratic variations, and second, 
our governing vector fields to include the associated second order partial dif¬ 
ferential operators shown above. Thus for i = 1 ,... ,d we set Di := Vi ■ d 
and 

x Vi] := [x l , X 1 ] and £> m := \Vi ® V: d 2 . 

Then the Ito chain rule takes the form 


f oY t = f oYq + Yj f D, 

aeA 


Of oYr dX° 


where A denotes the alphabet set of letters {1,..., d, [1,1],..., [d, d]}. Iterating 
this chain rule produces the formal Taylor series expansion for the solution and 
thus flowmap given by 

Vt = 'y \ lyj (f ) • 

w 

Here the sum is over all words/multi-indices w that can be constructed from 
the alphabet A. All the stochastic information is encoded in the multiple 
stochastic Ito integrals I w = I w (t ) while the geometric information is en¬ 
coded through the composition of partial differential operators D w . For a word 
w = ai ■ ■ ■ a n these terms are D w := D ai o ■ ■ ■ o D an and 

I w ■= [ dX%---dX%. 

JO^T 2 !S---^Tn^t 

It is natural to abstract the solution flowmap and view it as an object of the 
form 

^w®w, 

which lies in a tensor product of two algebras. The algebra on the left is asso¬ 
ciated with multiple integrals and the algebra on the right is associated with 
partial differential operators. The algebra on the left should be endowed with 
a quasi-shuffle product, to reflect the fact that the real product between two 
multiple Ito integrals I u and I v generates a sum over all multiple Ito integrals 
generated by the quasi-shuffle of the words u and v; we define this product 
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precisely, presently. The algebra on the right should be endowed with a con¬ 
catenation product, to reflect the fact that the composition of two differential 
operators D u and D v generates the differential operator equivalent to that 
represented by the concatenation of the words u and v. In the next section 
we define these underlying concatenation and quasi-shuffle algebras, their cor¬ 
responding Hopf algebras and the algebras associated with endomorphisms 
on them. These algebras prove useful in the following sections, they keep our 
proofs direct and succinct. 

3 Quasi-shuffle Hopf algebras and endomorphisms 

Our exposition here is based on Reutenauer [35], Hoffman m and Hoffman 
and Ihara [2(J. Let A denote a countable alphabet and KA the vector space 
with A as basis and K a field of characteristic zero. Suppose there is a commu¬ 
tative and associative product [ •, • ] on KA. We use K(A) to denote the non- 
commutative polynomial algebra over K generated by monomials (or words) 
we can construct from the alphabet A. We denote by A* the free monoid of 
words on A. 

Definition 2 (Bilinear form) We define the bilinear form (•, •): K(A) (x) 
K(A) —* K for any words n, v e A* to be 



This is equivalent to the scalar product given in Reutenauer [22 P- 17] and 
Hoffman p. 57]. For this scalar product, the free monoid A* forms an 
orthonormal basis. We will always assume that A equipped with [ •, • ] satisfies 
the following finiteness condition: for all letters c e A the cardinality of the set 
{a, b e A: ([a, b\, c) ¥= 0} is finite. It is satisfied, though not restricted to, when 
A is finite. The following example illustrates the case of a possibly infinite 
alphabet. 

Example 1 Consider a minimal family of general semimartingales given by 
{X 1 ,..., X d } in the sense outlined in Curry, Ebrahimi-Fard, Malham and 
Wiese [?]. We do not restrict ourselves here to continuous semimartingales. 
However, a collection of independent continuous semimartingales, or a collec¬ 
tion of independent Levy processes, is a minimal family. We can construct a 
countable alphabet A as outlined therein as follows. With each semimartin¬ 
gale we associate a letter 1 In addition, inductively for n ^ 2, we 

assign a distinct new letter for each nested quadratic covariation process 
[X k \[X k \[...[X k "-\X k "]...]]] with ki e for i = l,...,n, “pro¬ 

vided it is not in the linear span of {X 1 ,... ,X d } and previously constructed 
ones”. Due to commutativity the order of the letters in the n-fold nested 
bracket is irrelevant and associativity means that we can render all n-fold 
nested brackets to the canonical form of left to right bracketing shown, or more 
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conveniently [X fcl ,..., X fen ]. We denote the new distinct letters by [fci,..., k n \. 
Hence our underlying countable alphabet A consists of the letters 1,..., d and 
all possible nested brackets [ ■, • ] generated in this manner. For convenience 
we set [X 1 ] = X 1 and thus also [i] = * on K(A). 

We use K(A) to also denote the concatenation algebra of words with conca- 
tentation as product. If u and v are words in K(A), then their concatenation 
is uv e K(A). 

Definition 3 (Quasi-Shuffle product) For words u, v and letters a, b the 
quasi-shuffle product * on K(A) is generated recursively by the formulae: 
u*l = l*u = u, where ‘1’ represents the empty word, and 

ua*vb = (u* vb) a + (ua * v) b + (u * v) [a,b]. 

Endowed with this product K(A) is a commutative and associative algebra 
called the quasi-shuffle algebra which we denote by ^(A)*; see Hoffman |T9]. 
In the special case when the generator [ ■, ■ ] is identically zero on KA, it 
reverts to the shuffle algebra K(A) UJ of words with shuffle lu as product, 
where ua lu vb = (u lu vb) a + ( ua lu v) b. 

Example 2 The quasi-shuffle of the words 12 and 34 is given by 12*34 = 
1234 + 3412 +1342 + 3142 +1324+ 3124 + 1 [2, 3]4+ [1, 3]42 + 3[l, 4]2+ [1, 3]24 + 
13[2,4] + 31 [2, 4] + [1,3][2,4]. 

Example 3 A minimal family of semimartingales generates a quasi-shuffle 
algebra. This is proved in Curry et al. [7]. 

Definition 4 (Deconcatenation and de-quasi-shuffle coproducts) We 

define the deconcatenation coproduct A: K(A) —> K<(A) ® K<(A) for any word 
ui e K(A) by 

A(w) := ^ (uv, w)u0v. 

U,V 

We also define the de-quasi-shuffle coproduct A' : K(A) —* K(A) (g) K(A) for 
any word w 6 K{A) by 

A'(w) := y ^ I (u*v,w}u(g)v. 

IL,V 

The finiteness condition on A ensures that A' is well defined. Endowed with the 
concatenation product and de-quasi-shuffle coproduct K(A) is a Hopf algebra 
which we also denote by K(A). No confusion should arise from the context. In 
addition, when endowed with the quasi-shuffle product and deconcatenation 
coproduct K<A) is another Hopf algebra which we denote by K(A) # . The 
antipode in both cases is the signed reversing endomorphism. We denote by 
End(K(A) s|c ) the K-module of linear endomorphisms of ^(A)^. 

Definition 5 (Convolution products) Suppose X and Y are two linear en¬ 
domorphisms on the Hopf quasi-shuffle algebra K(A) # . We define their quasi¬ 
shuffle convolution product X * Y by the formula X * Y := quas o (X ® Y) o A, 
where ‘quas’ denotes the quasi-shuffle product on K(A) # . 
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Remark 1 We use the same notation for the quasi-shuffle convolution prod¬ 
uct as for the underlying product, no confusion should arise from the context. 
There is also a concatenation convolution product cone o (X g) F) o Z\' on 
End(K(A)), where ‘cone’ denotes the concatenation. 

In other words, since deconcantenation A splits any word w into the sum of 
all two-partitions u®v with u. v e A*, including when u or v are the empty 
word 1, we see that 


(X*Y){w)= ^ X(u)*Y(v). 

uv=w 

Now consider the following algebra which plays an essential role hereafter, 

K<A>*® K<A>. 

This is the complete tensor product of the Hopf quasi-shuffle algebra on the 
left and the Hopf concatenation algebra on the right; see Reutenauer 131 p. 18, 
29]. It is itself an associative Hopf algebra. The product of any two elements 
in this tensored Hopf algebra, which extends linearly, is naturally given by 

(u®u')(v ®v') = (u * u') (g) (vv'). 


Remark 2 In our context, the significance of this tensor algebra is that it is 
the natural abstract setting for the flowmap. 

Any endomorphism X e End(K(A) # ) can be completely described by the 
image in K(A) s .(g) K(A) of the map 

X ^ y X(w)®w , 

weA* 

see Reutenauer [31 p. 29]. Note that the identity endomorphism ‘id’ on K(A) 
maps onto J]w®w. Indeed the embedding End(K<A) H .) —» I^A)*® K<A) de¬ 
fined by this map is an algebra homomorphism for the quasi-shuffle convolution 
product. The unit endomorphism v on the algebra End(K(A) s .) sends non¬ 
empty words to zero and the empty word to itself. This embedding provides 
a mechanism for representing functions of J] w® w in End(K(A) ). Before we 
demonstrate this, we need the following. 

Definition 6 (Augmented ideal projector) We use J to denote the aug¬ 
mented ideal projector. This is the linear endomorphism on K(A) H . or K(A) 
that sends every non-empty word to itself and the empty word to zero. From 
the definition of the unit endomorphism v given above we see that J = id — 

We observe we can apply a power series function such as the logarithm function 
to the element ® w in K(A) # ® K<A) as follows. If T’ represents the 
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empty word, |w| represents the length of the word w and Ck '■= (—l) fc ^ for 
all k e N, then by direct computation we find 



w Q9 w 



k 



k 


2 Cfe 2 (mi * ••• *u k )®{u\---u k ) 

k^l m ,...,itfceA*\{l} 



(x) w 


where J* k denotes the fcth quasi-shuffle convolution power of the augmented 
ideal projector J. We emphasize the elements in the partition u\ ■ ■ ■ Uk = w in 
the sum in the penultimate line are all non-empty. Note that words cannot be 
deconcatenated further than all the letters it contains, and thus J*(w) is zero 
if w has length less than k. We conclude the action of the logarithm function 
power series on w(x)w can be represented by a corresponding power series 
endomorphism in End(K(A) # ). This will prove useful in Sections [4] and [5] so 
we summarize the result as follows. 

Lemma 1 (Logarithm convolution power series) The logarithm of the 
element w®w is given by 

log E ir(g)uu = ^ log* (id) o w (g) w, 

\we A* / u>eA* 

where 

f_i\k-1 

log* (id) := ^ J* k ■ 

l 

We will often abbreviate log*(id) o w to log*(to). 

We also note that equivalently, the embedding End(K(A)) —> K^A)*® K(A) 
given by 

Y i-» ^ w®Y(w), 

WE A* 

is an algebra homomorphism for the concantenation convolution product. 
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Definition 7 (Adjoint endomorphisms) Two endomorphisms A' and Y 
are adjoints if the images of 


X i—> ^ w ®X(w) and Y >—> ^ Y(w) 0 w, 


W W 


are equal. Hereafter we use X^ to denote the adjoint of X. 

This coincides with X and X”f being adjoints in the sense <(X'(zz),zA 


(u,X(v)y for all u,v e A*; see Reutenauer [Mj Section 1.5]. 

There is a natural isomorphism between the Hopf shuffle and quasi-shuffle 
algebras discovered by Hoffmann [191 which will play an important role here; 
also see Foissy, Patras and Thibon M for a theoretical perspective. To de¬ 
scribe the isomorphism succinctly we need to introduce the notion of compo¬ 
sition action on words. For any natural number n, we use C(n) to denote the 
set of compositions of n. A given composition A in C(n) will have say £ < n 
components so A = (Ai,..., Af). For such a A we define the following simple 
multi-index functions, ]A| := £ as well as 

17(A) := Ai + • • • T Xe, n (A) ■ Ai ■ ■ ■ A/ and T(A) ■ Ai! • • • A^!. 

We now define the composition action on words and the exponential map from 
Hoffman [19] . 

Definition 8 (Composition action) For a given word w = a± ■ ■ ■ a n and 

composition A = (Ai,..., A^) in C(n) we define the action A o w to be 

A o w := [oi • ■ • axJ[aA 1+ i • • • aAi+A 2 ] ■ • ■ [oAi+.-.+A^j+i 

where the brackets are concatenated. Here for any word w = ai ■ ■ ■ a n the 
notation [zc] denotes the n-fold nested bracket described above. 

Definition 9 (Hoffman exponential) We define the map exp H : IK(A) m — > 
K(A) # as that which leaves the empty word unchanged and for any non-empty 
word w is 



Its inverse log H : K(A) —> K(A) UJ for any word w is given by 



AeC(|iu|) y 1 


Hoffman m proved the exponential map is an isomorphism from the Hopf 
shuffle algebra K(A) UJ to the Hopf quasi-shuffle algebra K<(A) # . 

Example 4 In the case of the word w = 010203 the Hoffman exponential is 


given by exp H (aia 2 a 3 ) = aia 2 a 3 + |[a 1 , 02 ]a 3 + \ai[a 2 , a 3 ] + g[ai,a 2 ,a 3 ]. 
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The adjoint of the Hoffman exponential expjj: K(A) conc A , —» K(A) conc 5 is 
an isomorphism and defined explicitly as follows. Here we distinguish between 
two Hopf algebras with concatenation as product, we have K(A) A , with 
de-quasi-shuffle A' as coproduct and K(A) conc g with deshuffle S as coproduct. 
Then for any letter a from the alphabet A we have 

exp^(a) := 2 ^ ~ 

[ai,...,o n ]=a 

See Hoffman m for more details. Note expjj is a homomorphism for the 
concatenation product. Its inverse logy: K(A) conc s —* K(A) conc A , for any 
letter a from the alphabet as above is 

logjj(a) := ^ 2 - ai ... a,„. 

n^l [ai,...,a n ]=a 

We note K(A) LU A and K(A) conc s are dual Hopf algebras as well as K(A) # A 
and K<A> conCiZV .’ 


4 Exponential Lie series for continuous semimartingales 

We show the logarithm of the flowmap for a system of stochastic differential 
equations, driven by a set of d continuous semimartingales X( and governed 
by an associated set of non-commuting vector fields Vi, for i = 1 ,... ,d, can 
be expressed as a Lie series. We begin by emphasizing that for a given set of 
orthogonal continuous semimartingales, the only non-zero quadratic variations 
are those of the form [X 1 , X*] for i = 1 ,... ,d. In particular all third order 
[X*,X*,X*] and thus higher order variations are zero. Hence the generator 
[ •, ■ ], for example in the underlying quasi-shuffle algebra, is nilpotent of degree 
3. 

Remark 3 Examples of continuous (local) martingales are Brownian motion, 
time-changed Brownian motion and stochastic integrals of Brownian motion. 
The representation results of Doob, of Dambis and of Dubins and Schwarz, and 
of Knight (see Theorems 3.4.2, 3.4.6 and 3.4.13 in Karatzas and Shreve FTI] ) 
show that these examples are fundamental. Hence important examples for con¬ 
tinuous semimartingales are those just mentioned to which a continuous finite 
variation process, i.e. the difference of two real-valued continuous increasing 
processes, is added. 

We saw in Section[5]that the Ito flowmap has the form IwD w where the sum 
is over all words constructed from the alphabet A which contains the letters 
1,..., d as well as the letters [ 1 , 1 ],..., [d, d]. We proposed there, to represent 
the Ito flowmap by the abstract expression w ® w in K<A) # (x) K(A), the 
tensor algebra of the Hopf quasi-shuffle and concatenation algebras. We now 
make this more precise. Let I denote the algebra generated by multiple Ito 
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integrals I w = I w (t ) with respect to the semimartingales {X 1 ,... ,X d } or 
any nested quadratic variation processes generated from them, and constant 
random variable 1. 

Definition 10 (Ito word-to-integral map) We denote by /i: K(A) sl . —■> I 
the word-to-integral map /i: w >—> I w assigning each word we A* to the 
corresponding multiple Ito integral I w . 

The Ito word-to-integral map is a quasi-shuffle algebra homomorphism, i.e. we 
have n(u *v) = ) for any it, v e A*; see Curry, Ebrahimi-Fard, Malham 

and Wiese [7]. Let D denote the algebra of scalar linear partial differential 
operators that can be constructed by composition from the partial differential 
operators D*. 

Definition 11 (Ito word-to-partial differential operator map) We de¬ 
note by /2: K(A) —► IB) the letter-to-partial differential operator map p, : i i— * Di 
assigning each letter i e A to the corresponding operator Di. Recall the oper¬ 
ators Di are given, for i = 1 ,..., d, by Di := Vi ■ 8 and Don := ^Vi (x) V): 8 2 . 

The map /2 is a concatenation algebra homomorphism, i.e. we have p(uv) = 
fi(u)p(v) for any u, v e A*. Naturally p(g)/i: K(A)^®K(A) —> I®D is also an 
algebra homomorphism. With this homomorphism construction in place, we 
observe 

2 I W ®D W = (n®p)o ( ^ w ® w 

W ' W 

In principle we can compute the logarithm of u> ® w in K<A) # ® K(A) 
in the form of the quasi-shuffle logarithm as outlined in Section [3] However 
the basis in which we expand the flowmap and its logarithm corresponds, on 
the right hand side, to the terms D w in D which are compositions of the 
vector fields V) • 8 and second order partial differential operators ^Vi ® V): d 2 
for i = 1 The question is, how can we express the logarithm of the 

flowmap in Lie polynomials or in particular, in Lie brackets of vector fields? 
The natural resolution is to use the Fisk-Stratonovich representation of the 
flowmap. From the stochastic analysis perspective the procedure is as follows. 


Definition 12 (Fisk—Stratonovich integral) For continuous semimartin¬ 
gales H and Z, the Fisk-Stratonovich integral is defined as 

f Hr ° dZ T := f Hr dZ T + \[H, Z]t, 

Jo Jo 

where the ‘ ° ’ indicates Fisk-Stratonovich integration; see Protter |35) p. 216]. 

Lemma 2 (Ito to Fisk-Stratonovich conversion) Fori = and 

any function f: M. N —* K, for the integral term in the ltd chain rule we have 

fVi • d )f( Y r) dX r = \\ V i ■ d )f( Y r) " dX* - \[{Vi ■ 8)f(Y),X% 

JO Jo 
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and 


[( V;i • d)f{Y),X i ] t = [V, • 8) (Vi • 8)f(Y T )A[X i ,X i ] T . 

Jo 

Proof To establish the first result we set H = (V-i ■ 8) f(Y) and Z = X 1 in the 
definition for Fisk-Stratonovich integrals above. For the second result we use 
the Ito chain rule in Section[2]to substitute for (V) ■ d)f(Y t ) into the quadratic 
covariation bracket on the left. Then using that the bracket is nilpotent of 
degree 3 for continuous semimartingales, and zero if any argument is constant, 
establishes the result. □ 

Substituting the product rule 

(Vi ■ 8) (V ■ 8) = (V <g> V ): 8 2 + (V z ■ 8V) ■ d 

into the second result and that itself into the first result in Lemma [21 generates 
the Fisk-Stratonovich chain rule. 

Corollary 1 (Fisk-Stratonovich chain rule) For any function f : —> 

R we have 

f(Y) = f(Y 0 ) +2 [(V • d)f(Y r ) °dX* - ±2 (V • W) ■ 8)f(Y T )d[X\X%. 

i =i Jo i = i J o 

We emphasize, the differential operator in second term on the right is a vector 
field. In addition, the usual rules of calculus apply to Fisk-Stratonovich inte¬ 
grals. As in the Section [2] for the Ito case, we extend the driving continuous 
semimartingales and governing vector fields as follows. For i = 1,..., d we set 
X^A ;= [ X ’, X 1 ] as before, however we now set V yq := — \ (Vi ■ 8Vi) ■ 8. Then 
the chain rule above generates the Fisk-Stratonovich representation for the 
flowmap 

J'W V/J }, 
w 

where the sum is over all words w constructed from the alphabet given by A := 
{1,..., d, [1,1],..., [d, d]}. The multiple stochastic Fisk-Stratonovich integrals 
Jw = J w (t) are defined over the same simplex as for the multiple Ito integrals 
but with each nested integration interpreted in the Fisk-Stratonovich sense. 
The basis terms V w are compositions of vectors fields from the alphabet A with 
the assignment of each letter to each vector field as outlined above. We define 
the Fisk-Stratonovich word-to-integral map v. K^A)^ —* J by v. w *—> J w , in 
a similar manner to that for the Ito word-to-integral map. Here J denotes the 
algebra generated by multiple Fisk-Stratonovich integrals with respect to the 
semimartingales {X 1 ,..., X d ) or any quadratic variation processes [A 1 , A 1 ] 
generated from them, and constant random variable 1. Note the alphabet A* 
contains all the letters 1,..., d and [1,1],.... [d, d\. Let ¥ denote the set of 
partial differential operators constructed by composition from the vector fields 
V. and Von for * = 1,..., d. 
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Definition 13 (Word-to-vector field map) We denote by v\ K(A) —> V 

the letter-to-vector field map v\ i >-» V t assigning each letter i e A to the 
corresponding vector field V) and each letter [*, i] 6 A to the corresponding 
vector field Von , for i = 1,..., d. 

The Fisk-Stratonovich word-to-integral map v is a shuffle algebra homomor¬ 
phism while the word-to-vector field map v is a concatenation algebra ho¬ 
momorphism. Hence the natural abstract setting for the Fisk-Stratonovich 
representation of the flowmap is the complete tensor algebra K(A) LU (x)K(A). 
The map v ® v : K(A) ^tx) K(A) —> 1 ® V is an algebra homomorphism. 

The Fisk-Stratonovich and Ito representations for the flowmap must co¬ 
incide and so must their logarithms. The Fisk-Stratonovich integrals J w sat¬ 
isfy the usual rules of calculus, see Protter [33], while the basis terms V w 
are compositions of vectors fields. Hence the Chen-Stricliartz formula applies 
to the Fisk -Stratonovich representation for the flowmap, see Strichartz [36]. 
Recall from Lemma [T| the quasi-shuffle convolution logarithm of the identity 
log* (id) can be expressed as a power series in the augmented ideal projector J. 
When the quasi-shuffle convolution product reduces to the shuffle convolution 
product, also denoted lu, the corresponding shuffle convolution power series 
log m (id) is given by J — \ J m2 + ^ J m3 + • • •. 

Theorem 1 (Chen—Strichartz Lie series) The logarithm of the flowmap 
has the series representation 

lo gfe J wV w J = I] rfl J io S m (,w) V Mc 

where = v o log^w) and 

, uv <, v 1 -l ( \ +u (—l) d ( CT ) ^|cr| — 1 

log (w) = 2j c ° a i w ) mth c ° ■= —p— I 

creS|„| I* 7 ! V i i 



Here for any word w = a± ■ ■ ■ a n £ A* the basis terms are given by Vj-^u := 
[V ai , [V a2 , ..., [K„_i, V an ]h ■ ■ -]l]l, where [•, -] L is the Lie bracket, are Lie 
polynomials. The non-negative integers d(<r) denote the number of descents in 
the permutation a. 


Proof First, since the Fisk-Stratonovich multiple integrals satisfy the usual 
rules of calculus, the underlying product is the shuffle product lu; this is 
a special case of the quasi-shuffle product for which the quadratic variation 
generator [ ■, • ] is identically zero. Hence we can emulate the derivation of the 
quasi-shuffle convolution given in Section [3] to show that 


l°g Q] J w Vw) = log o(v ®h)o Q] w (x) w'j 
= {v (x) p) o log o f y w ® w^j 
= (v ® V) o log m (id) o w® u^j , 
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where the sums are over all words we A*. The expression log LU (id) is the 
shuffle convolution power series for the logarithm on the identity described 
above. That this power series has the equivalent expansion in terms of inverse 
permutations with the form for the coefficients c a shown, is proved in Section[5] 
for log*(id). See in particular Corollary^ Note that the coefficients «/ log uj( w ) 
then refer to the linear combination of multiple Fisk-Stratonoviclr integrals 
enumerated by the words generated by the set of permutations shown. 

Second, we express the logarithm of the Fisk-Stratonovich flowmap in 
terms of Lie polynomials as shown. The crucial observation here is that the 
adjoint of the shuffle convolution logarithm log m (id) is the concatenation con¬ 
volution logarithm on the identity (computed with the deshuffle coproduct 
<5) which we denote by log(id). This follows as for any integer k $= 1, the 
adjoint of J mfc is the fcth concatenation convolution power of J. The concate¬ 
nation logarithm of the identity is a Lie idempotent known as the Eulerian 
or Solomon idempotent. Another Lie idempotent is the Dynkin idempotent 
which has several characterizations, here it suffices to define it as follows. 
Let [1 • • -p] l denote left to right Lie bracketing of the word 1 • • -p so that 
[1 ■ • - p]l := [1, [2, [. . .[p - l,p]...]]] L . The element ^[1---p]l of K[§J is 
known as the Dynkin idempotent, where K[S P ] is the group algebra over K for 
the symmetric group S p of order p. We denote the Dynkin idempotent as 

6 := -[1 • • -p] L . 

V 

That the Eulerian and Dynkin operators are Lie idempotents is proved in 
Reutenauer [M* p. 195] for example. The Dynkin-Friedrichs-Speclit-Wever 
Theorem is now key. It states any polynomial P in the concatenation Hopf 
algebra K(A) lies in the corresponding free Lie algebra if and only if 9P = P. 
The image of the Eulerian idempotent is contained in the free Lie algebra 
associated with K<A>; see Reutenauer [3U p. 59]. Thus the element log(u>) for 
any word w e K(A) is a Lie element. The Dynkin-Friedrichs-Specht-Wever 
Theorem thus implies 

0 (log(u>)) = log(u>). 

Now by direct calculation, 

log (e W®Vj\ = »g-M 

= ^w<S) log(w) 

= ^w®6(log(w)) 

= ^ log m (u;) (x) 0(w) 

= ^n log m H <E> Ml, 

^ |w| 

where the sums are over all words we A*. Here we used that if X and Y 
are adjoint endomorphisms, and Z is an endomorphism on the concatenation 
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Hopf algebra, then we have 

^X(w) <S>Z(w) = Z w ® Z{Y{w)). 

We applied this identity with X = log m (id), Y = log (id) and Z = 0. □ 

The Hoffman exponential map exp H naturally relates Ito and Fisk-Stratonovich 
multiple integrals as follows; see Kloeden and Platen E31 Remark 5.2.8] for 
the Wiener process case. 

Proposition 1 (Ito to Fisk—Stratonovich: Hoffman exponential) For 

continuous semimartingales, for any word w in A*, we have J w = / e xp H 0) 
where explicitly we have 

T ex PhM = Z o^(A) —|A| t Xow = + Z ohFR^' 

AeC(H) ue[[iu]] 

The nilpotency of the generator [•, ■] implies the compositions X e C(|w|) with 
a nonzero contribution only contain letters 1 and 2 and so F(A) = 

Thus the set [[u>]] consists of the words we can construct from w by successively 
replacing any neighbouring pairs ii in w by \i,i\. 


Proof Using the definition of the Fisk-Stratonovich integral, for any word 
w = ai ■ ■ ■ a n we have 


Ja\---a n 


Jai—a„-i dd a „ + ^ d ai ... 0n _ 2 d[/ Qn _ 1 ,/ a „]. 


Recursively applying this formula for J Ql ... a „_ 1 and so forth, generates the 
result. □ 


With this in hand, we deduce the following main result of this section. 

Corollary 2 (Ito Lie series) Let c a denote the coefficient of a~ x {w) in the 
expression for log 1X1 (w) above. We can express the Chen-Strichartz Lie series 
in terms of multiple Ito integrals as follows 

l°g(Z J -M =Z Z 0^ ex PH( cr " 1 ( I "))^[' l "]L’ 

' w ' w <tgS| w | ' ' 


or equivalently, by resummation of the series, 


log ( Z JwVw ) = Z lw ( Z 


rUr 


[ct(w)]l 


Ao-eSi 


M 


+ z 

«e]Ml 


1 


2\ u \-\w\ 


z 


rVr 




treSi 


M 


Here for any wordw e A* the set ]]w[[ consists ofw and all words we construct 
from w by successively replacing any letter [*, i\ with ii, for i = 1,..., d. 
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Since the Fisk -Stratonovich and Ito representations for the flownrap must 
coincide, the expression above must coincide with the logarithm of the Ito 
representation for the flowmap. The key fact distinguishing the ltd from the 
Fisk-Stratonovich representation for the flowmap is that the word-to-vector 
field map v : K(A) conc s —> V and word-to-partial differential operator map 
p : K(A) conc A , —> ID), which are both concatenation homomorphisms, assign 


i ^Vi-d, ^ _ | i ~Vi-d, 

[M] \Vi ® Vi: d 2 , |[i, i\ >-*• -|(V) ■ dVf) ■ d. 


The ltd representation is constructed by composing the operators on the 
left shown above, while that for the Fisk -Stratonovich representation is con¬ 
structed by composing the operators on the right. The operators 5 V) ® V. : d 2 
and 2 (Vi ■ dVi) ■ d are both associated with the quadratic variation process 
[X 1 , X l ]. The ltd word-to-integral map ji : K(A) # A —»I and Fisk-Stratonovich 
word-to-integral map v. —» JI both assign i >->• X 1 and \i,i\ 

[X 1 , X 1 ]. The former is a quasi-shuffle homomorphism and the latter a shuffle 
homomorphism and consequently J w = I e xp H (w)- This relation, together with 
the calculus product rule given by (V^ • 5) (Vi ■ d) = (Vi® Vi): 8 2 + (Vi ■ dV) • d 
underlie the following result. Recall the definition of the Hoffman exponential 
and logarithm map adjoints in Section [3] 


Theorem 2 (Ito and Fisk—Stratonovich map relations) The two word- 
to-integral maps /i and v and the word-to-vector field and word-to-partial dif¬ 
ferential operator maps fi and v are related as follows 


v = fi o exp H 


and 


v = /i ° logn 


Proof The first relation follows directly from J w = I exPh ( w ) o v o w = 
p o exp H o w. The second relation follows from the calculus product rule above 
which can be expressed in the form 


\poii = fio[i, i]—vo[i, *] <»■ vo[i, i] = p,o ([*, i] — ^ii) o vo [i, i ] = ^ologjj o [i, i\. 


Note ii denotes the concatenation of i with i and p(ii) = v(ii). The final 
expression in this sequence follows using the nilpotency of the bracket [ •, • ]. 
Using that logjj o i = i for the letters i = 1,..., d and log^ is a concatenation 
homomorphism, establishes the second result. □ 
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Some immediate consequences of this result are as follows. First, algebraically, 
we observe 


2 I W D W = ^ j ^LOW® ^LOW 

= X ° w ° exp^ o w 

= (ji ® v) o w ® exp|j o w\ 

= (n ® v) o ^X ex Pn ow ® W J 

= X ^ ° ex Pn' 

= x^. 

In other words the Ito and Fisk-Stratonovich flowmaps coincide. Note the 
transfer of exp|j from the right of the tensor product to exp H on the left, relies 
solely on the vector space properties of K(A). Second, using this result we 
observe 


■ ow ®vo w 
v o w (x) V o w 


= [v®v) 


= [u®u 


,w®w 


ow®w 


log( X I wD w ) = log o(/z ® n) o l X - 

= log o{v ®v)o ( y, w ® vj\ 

') ° log O (X 

') ° (X lo s L 

= (v ® v) o ( V -L log m ow ® [ui]i 

M 

= (/r ® v) o ( V — exp H o log m ow ® [w]i 

M 

= X l^j^expnPog^^J^ML! 

where the sums are over all words w e A*. We can also perform a resummation 
of the series as indicated in Corollary [2] as follows, 

log (X = (u ® D) a (V log m ow ® 9 o w^j 

= (n ® v) o ( Y, exp H o log 11 " ow®9o w^j 

= (n ® v) o (X w®9 o log o expjj o«)j, 
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where 

log o expjj otc = log olw+ ^ ) 

V «6]]ll)[[ / 

= ^ c a a(w)+ £ ^ c a a(u). 

creSj ^ | we]]^[[ o-eS| u | 

These last two results are thus a restatement of the Ito Lie series results in 
Corollary [2] 

Remark 4 From the algebraic combinatorial computations above, we ob¬ 
serve: (1) In the first computation above the transformation from the Ito 
to Fisk-Stratonovich flowmaps was instigated by the transformation of coor¬ 
dinates fjL = v o exp^. In other words this transformation, which is a direct 
result of the product rule, encodes all the information required for Ito to Fisk- 
Stratonovich conversion; (2) Quadratic variations are a natural component in 
the Fisk-Stratonovich formulation; (3) The encoding which retains the let¬ 
ters [ 1 ,1 ],..., [d, d] as well as 1 ,..., d in the alphabet appears to be natural, 
especially in the context of using the quasi-shuffle machinery provided by Hoff¬ 
man m- Indeed this is also the case for stochastic differential equations driven 
by Wiener processes for which it is usual to replace the quadratic variation 
terms by the corresponding drift term; (4) The flowmap satisfies the linear 
equation ip = id + ^dS 1 with S := Yi i D i X' 1 , see Ebrahimi-Fard, Malham, 
Patras and Wiese HU. When the coefficients Di are constant, the solution is 
the well-known Doleans-Dade exponential; a representation of it in terms of 
iterated integrals was derived in Jamishidian m- 


5 Quasi-Shuffle Chen—Strichartz formula 

In this section we do not make any nilpotency assumptions on n-fold nested 
brackets as in Sectional We derive an explicit formula for the coefficients of the 
quasi-shuffle convolution logarithm of the identity endomophism on K(A) st , i.e. 
we explicitly enumerate 

log* (id). 

This represents the quasi-shuffle logarithm equivalent of the Chen-Strichartz 
shuffle logarithm formula and was derived in Novelli, Patras and Thibon [3D] 
and generalized to linear matrix valued systems in Ebrahimi-Fard, Malham, 
Patras and Wiese m- Using the notion of surjections instead of permu¬ 
tations and quasi-descents, we can closely follow the development given in 
Reutenauer [33j. We begin by outlining the theory of surjections and quadratic 
covariation permutations, which we hereafter call “quasi-permutations”, as 
well as their action on words. We denote the symmetric group of order p by S p 
and corresponding group algebra over the field K by K[S p ]. The crucial fact 
about any permutation a e K[S p ], which underlies the classical shuffle Chen- 
Strichartz formula, is that the inverse a records the following information: 
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“The letter i is at position cr _ 1 (z) in a”. We exploit the corresponding result 
for surjections herein. 

We shall denote the set of surjective maps from the set of natural numbers 
{1,... ,p} to the set of natural numbers {1,..., q} with q < p by We set 

% U %.r 

q^p 

Naturally we have S p cr §' Associated with each surjection in E>' p is a quasi¬ 
permutation. 

Definition 14 (Quasi-permutations) We denote by C(S P ) the set of all 
quasi-permutations, these are all the permutations in § p together with all 
unique words formed by applying all possible composition actions to these 
permutations, taking care to unify equivalent terms. In other words, 

C(S P ) := {A o a : <j e S p , A e C(p)}, 

where we identify all terms that are equal due to the symmetry and associative 
properties of the nested bracket operation. 

Henceforth we record quasi-permutations simply as a e C(S P ). However, it is 
always possible to decompose (non-uniquely) any given quadratic covariation 
permutation into its composition and permutation components, say as A o p 
or as a pair (A, p) where A is a composition and p a permutation. Given any 
quasi-permutation in C(S P ), there is a unique surjection in S p that records the 
position of letters and nested brackets of letters. 

Example 5 Consider the set of all quasi-permutations (£(§ 3 ), these are given 
by the set of §3 permutations 123, 132, 213, 231, 312, 321 and [1,2]3, 1 [2, 3], 
[1,3]2, 2[1,3], [2,3] 1, 3[1,2], [1,2,3]. The set of all surjections in §3 consists 
of 123, 132, 213, 312, 231, 321 and 112, 122, 121, 212, 211, 221, 111. Term by 
term in the order given, we see that the surjections record the corresponding 
positions of the letters in the quasi-permutations. 

Hence, by analogy with permutations, quasi-permutations play the role of gen¬ 
eralized permutations, while the corresponding surjections play the role of the 
inverse permutations by recording the positions of the letters in the corre¬ 
sponding quasi-permutations. Hence we have the corresponding statement to 
that above and crucial fact about surjections: each surjection ( corresponding 
to a given quasi-permutation a records the information: 

The letter i is at position £(z) in a. 


Example 6 The surjection 3221 from § 4 , tells us that the letter that was in 
position 1 in the quasi-pernrutation was sent to position position 3, the letters 
2 and 3 were sent to position 2, and the letter 4 was sent to position 1. Hence 
the corresponding quasi-permutation is 4[2,3] 1. For another example, if 2312 
is a given surjection in S 4 mapping 1 1 —»■ 2, 2 >—► 3, 3 >—»• 1 and 4 i—> 2, then the 
corresponding quasi-permutation is 3[1,4]2 which is equal to 3[4,1]2. 
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With each surjection we can associate a quasi-descent set. 

Definition 15 (Quasi-descent sets) Given any surjection ( e SJ, we define 
its quasi-descent set Des(£) to be the list of the indices k e {1,... ,p — 1} for 
which £(fc + 1 ) < £ (k). 

For the particular subset S p c S' p these inequalities would be strict and the 
indices k would correspond to the classical descent indices. Just as there is 
an intimate relation between shuffles and descents, there is also one between 
quasi-shuffles and quasi-descents. The following first key result underlies the 
whole of this section. 

Lemma 3 (Quasi-descents and quasi-shuffles) The set of surjections £ e 
S' p satisfying Des(£) c: {g} for q < p, is identical to the set of surjections 
satisfying £(1) < • • • < £( q ) and £(</ + 1) < • • • < £(p). 

Proof We observe that for any surjection £ e S' p for which Des(£) <= {g} for 
some natural number q < p, then discounting the case when the quasi-descent 
set is empty, by definition we must have £(fc) ^ £(fc + 1 ) ==> k = q k ¥= 
q =^> £(fc) < £(fc + 1). The latter condition is equivalent to that in the 

statement of the lemma. □ 

The second key result we establish in this section is a natural consequence. 

Corollary 3 (Quasi-descents and quasi-shuffles) Let q < p be natural 
numbers. If we factorize the word 1 • • -p = U\U 2 with |rti| = q and |zi 2 1 = p — q, 
then the quasi-shuffle product of u\*U 2 is given by 

u 1 *u 2 = Yj ct ( 0 “ Y f7 ^’ 

C(l)<— <C(q) Des(f)e{iir} 

C(9+l)<—<C(p) 

where <r(£) denotes the unique quasi-permutation associated with a given sur¬ 
jection £. The first sum is over all £ 6 S p satisfying the inequalities shown, 
and the second sum is over all £ e E> p such that Des(£) c {g}. 

Proof Recall the definition of the quasi-shuffle product and its generation 
through the formula 

■-q)*{q+l-■-p) = ((1 • ■ • — 1) * (<? + 1 • • ‘P)) q 

+ ({l---q)*{q+l---p-l))p 
+ ((!•■ -q -!)*(?+ 1---P- 1 )) [q,p]. 

We observe that if we recursively apply this formula to obtain on the right- 
hand side the complete sum over all quasi-permutations, then the quasi-shuffle 
product of zti and U 2 is equivalent to the prescription that it is the sum over 
all quasi-permutations whose corresponding surjections satisfy the set of in¬ 
equalities £(1) < • • • < £(g) and £(g + 1) < • • • < £(p). This establishes the 
first result. With this in hand, the result of the quasi-descent and quasi-shuffle 
conditions Lemma [3] above, implies the equivalence to the second result. □ 
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The following generalization is then immediate and represents the quasi-shuffle 
analog of Lemma 3.13 in Reutenauer [34] p. 65]. 

Corollary 4 (Multiple quasi-shuffles and quasi-descents) Let pi, 

Pk be positive integers of sum p and S = {pi,pi + P 2 , ■ ■ ■ ,Pi + • ■ • + Pk-i} 
be a subset of {1,... , p — 1}. If we factorize the word 1 • • -p = u± • ■ -Uk with 
|u,;| = Pi f or i — 1, • ■ ■) k, then we have 

ui*---*u k = 2 ^(C)- 

Des(C)SS 

Much like the symmetric group action on words there is an analogous quasi¬ 
permutation action on words. Recall we can decompose any quasi-permutation 
in K[C(S P )] as a = A o p, into its composition A e C and permutation p e 8 p 
components. 

Definition 16 (Quasi-permutation action) The action of K[C(S P )] on 
K(A) for any er e K[C(§ P )] decomposed as er = A o p and word w = aq ... a p , 
is defined by aw := A o (a p (i) • • • a p ( p )). 

We can now construct the quasi-shuffle logarithm of the identity. We start 
with the quasi-shuffle convolution powers of the augmented ideal projector J. 


Corollary 5 (Convolution powers and descents) For any k ^ 1 and 

word w we have 


J* k {w) = 2 E *(0 

V|S|=fc-l Des(C)s5 / 


o w. 


Proof For any word w the quantity J* k (w) is the sum over all possible k- 
partitions of w , say v\ ■ ■ ■ Vk, quasi-shuffled together. Hence we have 

J* k {w) = E Vi*...*Vk 

Vl ■■■Vk=W 

= E (ui * . . . *U k ) O w 


= E E -(C) 

|S| = fc-l \Des(C)SS / 

= ( E E -(o N 

V|S| = fe-l Des(C)ES / 

where we used Corollary |4] in the third step. 


o w 


o w, 


Corollary 6 (Quasi-Shuffle convolution logarithm on words) The ac¬ 
tion of the quasi-shuffle convolution logarithm on any word w is as follows 

i°g*M=( E Torrr E -(C)! ou ’- 

\S<={l J ... 1 |u)| —1} 1 1 " r Des(C)cS j 
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Proof By direct computation using Corollary [5] we find 

log» = £ 

1 



1)k -i / 

k \ 

£ 

£ 

k^l 


|S| = fc-l Des(C)ES 

(s 

(- 1 ) 1 *! 
|S| + 1 

£ 

o-(C) ) 


Des(C)£S 

/ 

\ss{r 

£ 

(- 1 ) |S| 

|S| + 1 

£ 

Des(f)E, 


o W 


O W 


o UK 


□ 


The following characterization of the quasi-shuffle convolution logarithm is the 
generalization of the standard shuffle convolution logarithm. We shall need the 
following integral identity for non-negative integers d and r which is proved 
for example in Reutenauer [3U p. 69]: 



(—1 ) d d\r\ 

(d + r + 1 )!' 


Corollary 7 (Quasi-Shuffle convolution logarithm endomorphism) 

The quasi-shuffle convolution logarithm log* (id) acts on 1 • • -p as follows 


iog*(i •■•/>)= £ 


(-rffio 

p 



o-(C), 


where d(£) denotes the number of quasi-descents in 


Proof We observe from Corollary [ 6 ] that log*(l • • ■ p) consists of a linear com¬ 
bination of quasi-permutations <r ((). Hence we directly compute the coefficient 
of an arbitrary quasi-permutation cr(C) in log*(l • • -p) which, using the result 
of Corollary [ 6 l is given by 


£ 

5e{l,...,p-l}: Des(C)£S 


(- 1 ) |S| 

|5| + f 


Suppose that ( has quasi-descent indicies pi,...,pk so that d(C) = k. To 
compute this coefficient we therefore have to determine the number of subsets 
S c {1,.. ,,p — 1} which contain pi,...,pi-. Note the coefficient itself only 
depends on the size of such sets. These subsets have possible size |S'| = k 
through to |S'| = p — 1. Starting with the case |S| = k there is of course only 
one set of this size containing pi,... ,pk, the set of these integers themselves. 
Now consider the case |S| = k + 1. Then an extra “quasi-descent” can be 
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placed in total of p — 1 — k possible positions, or equivalently in p — 1 — k 
choose 1 ways. When |£| = k + 2, there are p — 1 — k choose 2 ways, and so 
forth so that in general, when IS] = k + i, there are p — 1 — i choose i possible 
ways. Hence the coefficient above equals 

-l-d(CA (~l) d(0+i _ 
i ) d(£) + * + 1 ' 

This form of the coefficient of is equal to 

J° cc d(c) (l+a;) p - 1 - d(c) dx 
(_i)d(C) d (C)!(p — 1 -d(Q)! 

p\ 

(_l) d (C) , p _ A - 1 

p U(o ) ’ 

using the integral identity preceding the corollary. □ 

Remark 5 This is equivalent to the quasi-shuffle logarithm given in Ebrahimi 
Fard et al. m Theorem 6.2]. We included a self-contained derivation here for 
completeness. 



p-l-d(C) 


E 

1=0 


6 Concluding remarks 

There has been a recent surge in the development of quasi-shuffle algebras 
and stochastic Taylor solution formulae in the context of semimartingales, on 
the theoretical and practical level. See Platen and Bruti Liberati [32], and 
Marcus [25], Friz and Shekhar m and Hairer and Kelly [T 8 | for contempo¬ 
rary references. For example Li and Liu [23; considered systems driven by 
both Wiener and Poisson processes. We have shown that the Chen-Strichartz 
flowmap solution formula which is well-known for Stratonovich stochastic dif¬ 
ferential systems driven by Wiener processes extends to systems driven by 
general continuous semimartigales. We demonstrated it is in fact a Lie series, 
and this property holds irrespective of whether we consider the system in the 
Ito or Stratonovich sense. We also give and prove an explicit formula for the 
Lie series coefficients. Curry, Ebrahimi-Fard, Malham and Wiese [7] have de¬ 
veloped so-called efficient simulation schemes for such systems driven by Levy 
processes. This involves the antisymmetric sign reverse endomorphism rather 
than the quasi-shuffle logarithm endomorphism. 
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